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Abstract 

A small difference between the energy of the para-positronium (p-Ps) and ortho-positronium 
(o-Ps) states suggests the possibility of superposition of p-Ps and o-Ps during the formation of 
positronium (Ps) from pre-Ps terminating its migration in the matter in a void ('free volume'). 
It is shown that such a superposition decohere in the basis of p-Ps and o-Ps and the decoherence 
time is estimated. The time scale of the decoherence estimated here motivates respective correc- 
tion in decomposition of the positron annihilation lifetime spectra. The way of the correction is 
sketched. The timescale of the decoherence suggests a need of awareness when experimental data 
from positron annihilation techniques are processed. More generally, the superposited state of Ps 
should contribute to the evolution theory of positronium in matter. 
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I. INTRODUCTION 



The most powerful model describing the formation of the positronium atom in the mat 



1]. A positron from the irradiation source 



2] 



is less than 



ter was developed by Stepanov and Byakov 

(the kinetic energy of positron produced by the most exploited isotope 22 Na 
0.5 MeV) passing through the sample creates a set of electrons, ions, radicals etc. losing 
its own energy. According to this model, at the beginning, each of the electrons caused by 
ionization receives about 1 keV of the kinetic energy and is able to ionize other particles. 
Further ionizations release 30-100 eV each and produce new charged particles which form 
clusters called spurs. If the energy of the positron decreases to about 1 keV, the spacings of 
ionization spurs become less distant and overlap each other. Finally, the positron ends up in 
the spreading cloud of ionized particles called blob with the energy lower than the ionization 
level. Further loss of the positron energy is transformed into vibrations of molecules and 
reaching thermalization energy the production of positronium atom arises as a profitable 
process. 

Before the formation of positronium atom, a positron and an electron can leave the blob as a 
weakly interacting, delocalized system (pre-Ps) migrating to near free volume where it forms 
the localized Ps. It is usually believed that at this time the whole population of Ps consists 
of p-Ps (spin S—0) and o-Ps (S—l) so the spin value for each particle is definite. Such a Ps 
can annihilate in different ways: 1). the annihilation of p-Ps produces two gamma quanta, 
2). o-Ps produces three gamma quanta, 3). each Ps can annihilate with an electron captured 
from the environment [pick-off process). The probability of pick-off annihilation increases 
with the local electron density and dominates in the case of annihilation in condensed mat- 
ter. In the case of o-Ps the annihilation via pick-off is a two-quantum process. However, 
when the free volume is not reached, the positron annihilates as a quasi-free particle with 
one of the electrons from the bulk (in particular those from the blob). 
The most serviceable channel of the positron annihilation in a matter is long-living o-Ps 
which succeeded in reaching the free volume and annihilating by pick-off. In order to sepa- 
rate the long living component from the lifetime spectra, they have to be decomposed into 
elementary channels of annihilation (free positron annihilation, ortho- and para-Ps annihi- 
lation, etc.). The exponential form of annihilation for each channel is commonly assumed. 



II. MOTIVATION 



The papers exploiting the blob model assume the definite value of spin number for pre-Ps 
(see positronium formation formulas in However, as the difference of energies of p-Ps 
and o-Ps states in the vacuum is about 8.4- 10~ 4 eV one can expect that the spin quantity 
could change during the migration of a pre-Ps through the sample due to the intermolecular 
interactions. It seems to be an oversimplification considering the definite spin value only, 
e.g. only o-pre-Ps. Instead, one should suppose that quasi-free Ps enters a free volume as 
a superposition of para and ortho states and thus the spin number could not be a proper 
number for such initial state. 

The value of Ps spin plays an important role in further analysis of Ps life. The annihilation 
of Ps, as an electromagnetic process, should obey the principle of conservation of the parity 
(here charge parity). This number for Ps in its ground state with the spin value S is 
7T(7 = (— l) s whereas for n photons it is ttq = (— l) n [5]. So until the spin number is 
well defined, the system cannot decay with the emission of any definite n photons. Thus 
supposing that Ps is the superposition of p-Ps and o-Ps after entering free volume, one 
could expect that Ps can "live" in the free volume without giving an identifying decay signal 
(the number of photons is undefined) until it becomes p-Ps or o-Ps. This implies that one 
cannot decompose the PALS spectra simply into the exponential curves for p-Ps and o-Ps 
because the decay (two-quantum, for example) cannot start just after the birth of Ps. In 
our opinion, as long as there exists the population of Ps which is superposition of spins, 
one should make the correction when decomposing the PALS spectra (the analysis of the 
correction way is not the subject of this paper but an example of applying the correction 
is shown in Appendix |A"|) . In this paper we show that it is expected to decohere Ps in the 
p-Ps and o-Ps basis and we estimate the time of decoherence. 



III. RATE OF DECOHERENCE. METHOD OF CALCULATION 

The most convenient hamiltonian to consider the spin system interacting with the en- 
vironment of spins (bath) is the Heisenberg hamiltonian p] which is the phenomenological 
expression incorporating an overlapping of wave functions of electrons of the system and an 
environment [?], S]- The problem of the decoherence of certain spin systems, in particular 
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consisting of a pair of particles interacting with the bath, was considered in 9KL3J which 
support our calculation method. 

Let us consider the hamiltonian for interaction of Ps and the molecules from the wall of the 
free volume 

2 N+2 
i=3 

where J- coupling constant, a^- Pauli matrices for the electron in Ps, erW- Pauli matrices 
for the i-th (of N) electron from an environment, the sign o denotes a scalar product. 
We assume the existence of free volumes where Ps could live sufficiently long timej. As 
the environment we mean atoms which form the wall of the free volume. There is a wide 
spectrum of materials where such free volumes exist (molecular crystals or polymers). In all 
these cases we suppose the atoms have certain electron density around (inside) the free vol- 
ume. These valence electrons can interact with an electron from Ps by exchange interaction 
and resulting in the pick-off process of Ps. 

The radii of the free volumes are in the range considered by the PALS method, i.e. from 
about lA to hundreds of nanometers. The wavelength of valence electrons lies in the range 
of angstroms. So one can expect in general that Ps in the free volume interacts with the 
electrons only momentarily. Generally, it is not obvious how many (N variable in this work) 
electrons interact with Ps at the same time. 

Suppose the initial (just after the localization of Ps in the free volume) state of the whole 
system is p = l^o) (^o|? where 

N+2 

l*o> = \Ps) J] \ Si ) = \Ps) \S) (2) 

i=3 

and 

\Ps)o = |(|0,0> + |1,-1) + |1,0) + |1,+1». (3) 

The ket \S) denotes the state of electrons of the environment whereas the numbers in the 
kets are spin s and its z-th projection s z for Ps, respectively. The latter equation expresses 
the supposition that after many interactions with the bulk, Ps, just after entering into the 
free volume, does not distinguish any particular Ps state. 



1 Our model seems to be the most suitable for the solid state of matter. Probably the model could be 
applied also to the liquid state where, according to the recent knowledge, the bubbles are formed as free 
volumes for Ps. 
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The evolution in time gives the following expression for the whole state at the instant t 

p(t) = \%)(9 t \, (4) 
where \^t) = e~ lHt ' n \^o). For the positronium only, we calculate the state at t as 

p Ps {t) = Tr env [p(t)} (5) 
tracing over the environmental space. 

In order to estimate the admixture of p-Ps in o-Ps at t we calculate some nondiagonal 
elements of pp s (t). More precisely, we consider the expression 

S(t)= jr \(s = 0\p Ps (t)\s = l,s z )\ 2 (6) 

s*=-l 

which is expected to collapse with time. The time when it reaches the minimum value we 
define as the time of the decoherence, 

Because of the mathematical problems with calculating the whole p(t) analytically, we made 
the following simplification 

|* t ) = e- iH ^\%) ~ V-^-) 3 l*o>, (7) 

j=0 3 ' 

where the value n is fixed in this formula. This truncation at appropriate n still allows to 
approximate T d sufficiently for our aims (see Appendix [Bl. 



IV. RESULTS AND DISCUSSION 

The numerical calculations of ([6]) using (|7j) bring to the following statements: 

1. Indeed, one can observe the decrease of S(t) in the basis of p-Ps and o-Ps states (fig. [Tj), 
but nondiagonal terms do not disappear completely and permanently. 

2. After the decrease of S(t) shown in figs [T] and [2] the oscillations were obtained^) which are 



consistent wit 



spin bath 10, 



i the results of calculations of decoherence for other complex systems in the 
11]. For N=2, Td ijf| about 29 ps which means that the effect of the existence 



2 Oscillations are visible when calculated with sufficiently large n to obtain reliable outcomes for t > Td. 
Here maximally n=10 was applied which is enough to calculate reliably the evolution up to about 30 ps, 
so reliable oscillations are shown for the settings with enough small Td- As an example see the calculation 

for iV=4, n=10 shown in fig. [2] and in figs |4] and [5] 

3 This result is not shown in any figure. Fig. [T] shows only a part of the whole S(t) for this N in order to 

show that each nondiagonal term constituting S(t) collapses jointly but differently. The collapse of the 
whole S(t) is shown in the next figures. 
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FIG. 1. |/0[o,O;i,+i]W| 2 calculated for a random initial state vector ([2]). The interaction of Ps with 
N=2 electrons from the wall was assumed and n in ([7]) was set 10. The time to the minimal value 
of the function is regarded as the decoherence time. 
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FIG. 2. S(t) averaged over 100 randomly chosen samples of the initial state vector for different 
electrons interacting with Ps atom (N). The comparison of calculations with truncation of the 
evolution operator by n=3 and by n=10. The results are shown for N=l (solid and dashed line, 
respectively) and for N=A (dotted and dashdotted line, respectively). 




of the superposition of states lasts for the time detectable in PALS experiments and one 
should take into account this effect when decomposing the PALS spectra. This statement 
we consider as the most significant conclusion of this work. 

3. Higher orders of hamiltonian in the evolution operator do not affect substantially the 
value of the decoherence time, within the meaning that Td may vary by few picoseconds, 
i.e. the decoherence maintains always the same time scale (fig. |2]). Figure |2] compares the 
calculations with the 3-rd and the 10-th power of H as the highest order for a given N. 
Skipped parts of expansion of the evolution operator may change Td to some extent but 
their contribution is estimated in Appendix [Bj The aim of this paper is to estimate Td but 
not to calculate the time evolution of the system at any time. The calculation is performed 
with sufficient accuracy required to estimate the values of the Td for a given N . To cal- 
culate further evolution, next powers n are needed. For example, in order to extinguish 
the increasing value of polynomials for the increasing time, t > Td, one needs to calculate 
next powers of H in the expansion of the evolution operator which give higher powers of 
polynomials dominating at large values of t; thus one obtains oscillations at times larger 
than Td instead of a rapid increase. Because the omitted terms of expansion may shift S(t) 
for larger values of time, these terms may be more appropriate correction of Td for slow 
decoherence (small number of electrons in the bath giving larger Td). 

4. Decoherence time depends on a coupling constant J, which was set assuming that the 
interaction energy of Ps and the electrons from the wall is about 10 -4 eV (this is the energy 
which could convert ortho and para states and cause the initial superposition state in the 
bulk. Such a value of energy for the interaction of Ps in the free volume is reasonable 
only for some specific cases, e.g. when ortho-para conversion takes place. This process is 
characteristic only for specific substances {l2j]). If the energy of the interaction is smaller, 
one can expect that the decoherence time increases (e.g. if the interaction energy is set to 
10~ 5 eV then Td ~150 ps for N=3 instead of Td=18 ps as estimated above). 

5. The value of Td decreases with the number iV of electrons in the environment reaching 
the asymptotic value 12.2 ps (fig. ED - 

6. If one assumes in the numerical calculations the magnetization of the medium (all spins 
in the same direction), like paramagnetics at extremally low temperatures, one obtains the 
suppression of the decoherence: it makes Td shorter for a given N and the minimum of 
S(t) becomes shallower, the shallower the greater is N. Figure H] shows the comparison of 
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FIG. 3. The decoherence time as a function of the number N of interacting electrons with Ps. 
The value is averaged over 20 randomly chosen initial state vectors for n=7. 




5 10 15 20 

tips] 



FIG. 4. The presence of the magnetization of the medium. The curves are calculated for N=4, 
n=10. The curve without magnetization (randomly situated spins of the environment) was obtained 
as an average over 100 random samples of spin set; dashed line. The bold line- all spins are situated 
along the z-axis, the dotted line- the direction of all environmental spins are situated along the 
x-axis. 
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FIG. 5. The case of magnetization along the z-axis. S(t) for N=l, 2, 3, 4 (lines: regular, dashed, 
dotted, dotdashed, respectively); n=10. 

an averaged S(t) for random initial vectors of the state (lack of magnetization) and for all 
spins 'up' in the initial vector (presence of magnetization). In the case of magnetization the 
minimum is shalloweiQ. Furthermore, fig. [5] shows S(t) for different iV for magnetization. 
One can see that for greater N one obtains shallower minima. 

Calculations concerning magnetization were performed when all spins are situated along 
different direction from that we considered Ps state. The result is shown in fig. HJ A sys- 
tematic study of Td dependence on the direction of electron spins was not performed. 
The effect of inhibiting the decoherence in the presence of magnetization along the z-axis 
is for us a new one and we regard it as a hypothesis to be proved. The situation resembles 
the Ps eigenstates in the magnetic field - it is known that in the case of external magnetic 
fields the eigenstates of respective hamiltonian are those which are mixtures of p-Ps and 
o-Ps. Here supposing the magnetization of the sample one introduces a magnetic field. For 
this case the superposition of p-Ps and o-Ps seems to be maintained (which is indicated by 
the shallow minimum). 

The case of the magnetization presented above cannot be treated as one for the ferromag- 

4 We decided to compare in this figure an averaged vector representing the state without magnetization 

with the state with magnetization where all spins of the environment are 'up'. The reason for averaging is 

that particular states of the environment differentiate the state of the whole system to some considerable 

extent and the only way to compare the results with magnetization is to average them. 
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netic sample. In ferromagnetic crystals Ps does not form. Such property of the bath may 
be rather connected with those of some molecules with the magnetic moment. 
7. The Heisenberg hamiltonian leads to reasonable results in the case of positronium deco- 
herence problem but it shows that the decoherence is not complete in the basis of p-Ps and 
o-Ps (there exists at any time the population of Ps atoms in the superposition state). The 
residue means that there exists at this time a Ps population with not well defined spin value 
or, in other words, not all Ps atoms are in the spin state. The probability to find such an 
undefined spin can be calculated as 1 — (\{p — Ps\Ps(t))\ 2 + \{o — Ps\Ps(t))\ 2 ) (it is neither 
p-Ps nor o-Ps). The smaller is S(t) value at minimum, the more Ps have finally well defined 
spin. 

We believe that the existence of the residue may be an artefact of the Heisenberg hamil- 
tonian being the simplest way of treatment of the realistic case. We intend to follow the 
calculations with a more complicated hamiltonian, where J is time dependent (Ps atom 
does approach to or goes away from the wall, so the interaction strength changes). 



V. CONCLUSIONS 

Our theoretical study of the decoherence of Ps in matter results in statement that co- 
hered Ps transforms into the states with well defined spin and that the decoherence lasts 
enough time (picoseconds and more) to modify the positron lifetime spectra and momentum 
distributions and should be taken into account when any interpretation of the results in 
positron studies is performed. We have shown the way of modification for formulas used for 
PALS spectra decomposition. 

The time of the decoherence decrease exponentially with the number of electrons interacting 
with Ps in the free volume of matter. We examined the influence of magnetization to be 
expected on the process of the decoherence. 

Our calculations are consistent, according to the formalism, with the previous ones pub- 
lished for similar quantum systems and the results add the new aspect to the recent theory 
of Ps formation and its evolution in matter. 
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Appendix A: 



The influence of hypothetical coherence of Ps on the PALS spectrum seems to be a 
complex problem which requires more study. As an example where this influence can be 
observed we propose the following one. 

To interpret the PALS spectra one considers the formula for coincidence (2- and 3- quantum) 
rate (iV 7 ) resulting from the formula (see also [3]]) 



dA 7 ~ -dP e (Ae) - dS { Q K) - dS { p ^ ] - d{ P Ps) {x " Ps) - d{oPs) {x ° p ° 



(Al) 



where: P e - amount of free positrons, S p , S Q - the same for pre-para-Ps and pre-ortho-Ps, pPs, 
oPs- amount of para- and ortho-Ps, respectively; As are the decay rates for a particular type 
of particles. 

Furthermore, one assumes P e , S p , etc. obey the following equations 



dP e {t) = -(A e + v)P e (t)dt, P e (0) = 1, 

dS (t) = \vP e {t)dt - (A e + K)S Q (t)dt, S o (0) = 0, 

dS p (t) = \vP e {t)dt - (A e + K)S p {t)dt, S p (0) = 0, 

d(oPs)(t) = K S {t)dt - X oPs (oPs){t)dt, (oPs)(0) = 0, 

d(pPs)(t) = K S p (t)dt - X pPs (pPs)(t)dt, (pPs)(0) = 0, 



(A2) 



where u, A e , etc. are the constants describing creation and decay rates of each population. 
If one considers an annihilation in matter where the pick-off process dominates, flAll) relates 
mainly to the 2-quantum process of annihilation and both flAlj) and flA2j) allow to calculate 
the intensity of the annihilation rate for the set of particular channels (i.e. for P e , (oPs), 
etc.). The decomposition of PALS spectra is a direct consequence of the use of formulas 
above. 

If one supposes the existence of the coherence of the states, one does not distinguish the spin 
value for pre-Ps and one assumes the existence of Ps (and also pre-Ps) in the superposition 
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state, which annihilates with the rate Am- So, one can replace (1A2|) and (lAlj) by 
dP e (t) = -(A e + v)P e (t)dt, P e (0) = 1, 

dS(t) = uP e (t)dt - (A e + K)S(t)dt, S(0) = 0, 
d(oPs)(t) = K S(t)dt - X oPs (oPs)(t)dt, (oPs)(0) = 0, 
d{pPs)(t) = K S(t)dt - \ pPs (pPs)(t)dt, (pPs)(0) = 



(A3) 



and 



dA" 7 ~ -dPj: Xe) - dS {Xe) - d(MS) {XMs) - d{pPs)^ Ps) - d(oPs) {x ° Ps) - d(MPs] 



i^MPs) 



(A4) 



The formula (IA4j) is not a simple replacement of (lAlj) because the last term in it describes 
the decays where the number of photons per one annihilation act is not fixed and, for 
example, 3 quantum decay cannot be neglected here. 



Appendix B: 

The H n operator can be written as 

N+2 n-fold sum 

H n = ( J2a®oaM) n = ^oa^-a^oa^-.... (Bl) 

p=3 p,s,t,... 

Each superscript changes in the range (3, N + 2), so we have N n different terms with the 
fixed value of indices. Because 

— J2^=i a P " cr | P) we nave 3" different values of each 
term summed in (IBip (for fixed p, s, t, . . .). Each of such terms can give maximally the 
value 1 on any stataj so the resultant coefficient related to the n-th power of hamiltonian in 
([7]) is K/n! • (—iJh/4 ■ t) n , where K = (3 AQ n . Substituting the assumed value of J into this 
coefficient we can see the contribution of the n-th term as a function of t. This contribution 
for fixed £=10 ps and t=20 ps (the region of for several N) is shown in figure |6j 
If one knows the contribution of skipped terms in the evolution operator at any t one can 
estimate the expected shift of minimum (T^) of polynomial we obtained from truncation at 
smaller n and that one expected when more accurate calculations could be performed. For 
doing so, see Appendix O 

5 We assumed a very coarse approximation: if one of gave 1 in the extreme case, then none of for 
i j could reach this value. In spite of this such extreme value is applied for simplicity. Because of this 
simplification the approximation is overestimated. 
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t=10ps 



t=20 ps 




FIG. 6. The boundary value of the absolute of n-th term (normalized to the maximal value) in ([7]) 
showing its contribution at t=W ps and i=20 ps to approximate the exact result in the region of 
calculated T^. 

Appendix C: 

Let us suppose we have two polynomials fi(t) = Yli=o a * ^ an< ^ fz(t) = Si=o k t l which 
represent the S(t) function calculated when n=k and n=k J r c respectively. Let single out 
s-th minimum of each polynomial, t ls , t 2s , and relate it to those constituting T d . Suppose 
we know (see Appendix [B]) amount of correction at any t given by c additional terms of 
expansion of evolution operator; denote this surplus as 5(t). This function S(t) can be 
expressed as 

k k+c 

S(t) = f 2 (t) - h(t) = - a t )f + ( C1 ) 

i=0 i=k+l 

Let us construct the auxiliary function 

F = f 2 (t)-f l (t-At), (C2) 

which has, if t=t 2s , a maximum for At=t 2s — t± s . Inserting the definitions of fi(t) into (1C2|) 
and using the Newton formula for the binomial expansion we have 

k 

F = 6(t) + Y}iai f- 1 At - ai l -^^t l -\At) 2 + ■ • ■ ] (C3) 

i=0 

and finally inserting into (IC30 the value t 2s we find a local maximum of F as 

k 

S(t 2s ) + Y}i ai t^At - a, l S*-A? 2 -\At) 2 + ••■]= max. (C4) 

i=0 
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This equation is equivalent to the derivative of its left side over At equated to 0. This 
equation should give a maximum value for At = t 2s — ti s , i.e.: 

dF 

7^|{t=t 2s , At=t 2s -tis} — 0- (C5) 

The polynomial f\ (t) is usually known as it is the solution of our problem for truncation of 
the evolution operator at smaller n. For known values of t± s and one obtains an algebraic 
equation for t 2s . In this way we can calculate the expected difference between known T d , 
obtained for shorter expansion of the evolution operator, and that one if more accurate 
calculations could be performed. 
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